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Abstract 

We show that the master equation governing the dynamics of simple diffu- 
sion and certain chemical reaction processes in one dimension give time evolu- 
tion operators (Hamiltonians) which are realizations of Hecke algebras. In the 
case of simple diffusion one obtains, after similarity transformations, reducible 
hermitian representations while in the other cases they are non-hermitian and 
correspond to supersymmetric quotients of Hecke algebras. 
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It is well known in the literature that several integrable quantum chains corresponding 
to magnetic systems U can be represented as generators of Hecke algebras H n (q), 
when certain artificial interactions are added in the bulk and the surface. In this letter 
we show that the master equation describing the dynamics of some chemical processes 
limited by diffusion give representations of Hecke algebras, where the supplementary 
interactions appear naturally. 

The Hecke algebra H n (q) (with n = L — l) is an associative algebra with generators 
ti (i — 1, . . . , L — 1) satisfying the relations 

eie i±1 ei - e, = e i±1 eje i±1 - e i±1 (1) 
[e i>ei ] = ; \i-j\>2 (2) 

^ = (q + q- 1 ) e, (3) 

where q is in general a complex parameter. In our applications q is always real. The 
main interest in this algebra in statistical mechanics is due to the fact that a spectral- 
dependent R(u) matrix, satisfying the Yang-Baxter relations 0, can be constructed 
in a standard form. This construction, also called "Baxterization" H], implies that 
the Hamiltonian J2f=i e % has an infinite number of conserved charges and this is the 
first step towards the exact integrability of H . 

One can define [3| a sequence (P, M) of quotients algebras of H q {n) (P and M are 
positive integers). Since the difinition is quite complicated, we will give only a few 
examples. The quotient (1, 0) implies taking e\ = 0. This quotient is not relevant for 
our purposes. The (2, 0) quotient is the Temperley-Lieb algebra: 

eie i± iei = e { . (4) 

The (1,1) quotient reads: 

(eie i+2 ) e i+ i (q + q' 1 - e^j (g + q' 1 - e i+2 ) = 0. (5) 

For our purposes it is sufficient to mention that each (P, M) quotient has as a repre- 
sentative one of the (P, M) Perk-Schultz || quantum chains where the Hamiltonian 
is given by the sum of H n (q) generators: 

H {P,M) = Y^ e {P,M) (6) 
3=1 



where 

(P,M) _ 'I - <l I v 



e 



i + 1 ' ( p«/3p^ , q + q 1 y^ 1 c ™ r, 

3 - o Z^i + o 2^ €ah 3 L 3+l 



o=0 



-,-1 



+ s —z— e ««« - m aE ?^ ] • (7) 
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with N = P + M. The N x N matrices E a/3 have elements: 



{E°*) yS = S^ps (a,fl = 0,l,...,N-l) 
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and 



6l — ... — Ep-\ — —€p — ... — — €p + M-l- 



(9) 



The quantum chains if( p > M ) are U q (SU(P/M)) invariant @ and have the important 
property that they contain all the irreducible representations of the (P, M) quotient. 
In the examples we are going to mention, if a chain corresponds to a certain quotient 
(P, M) them its energy levels (not the degeneracies!) are the same as those of the 
representative of this quotient given by Eqs.(^|-0). We now turn to the physical 
problem. 

It was pointed out recently |7[], that the one-dimensional master equation corre- 
sponding to a chain with two-body transition rates can be written as a Schrodinger 
equation with a Hamiltonian described by a quantum chain with nearest-neighbour 
interactions. Several examples were found in which this Hamiltonian had precisely 
the structure of Eqs.(^-|7]). In this letter we present this problem in detail. We will 
use the notations of Ref. 0. 

In order to describe the dynamics of general react ion- diffusion [7j processes we 
use a master equation. Let us consider an open chain with L sites, where at each site 
i — 1, 2, . . . , L we attach a variable (3 taking N integer values (/3 = 0, 1, . . . , N—l). By 
convention we attach the value (3 = to a vacancy (inert state) and (3 — 1, . . . , N — 1 
a molecule of type (3. The master equation describing the time evolution of the 
probability distribution P( {/?};£) is 



where we, m are the transition rates and the prime in the second sum indicates the 
exclusion of the pair I = m = 0. In ([IT]) the symbol [x + y]jv means the addition 
(x + y), modulo N. We restrict ourselves to the cases where the transition rates are 
homogeneous and depend only on the molecules on nearest neighbouring sites. The 
transition rate w^ m (a, 0) gives the probability, per unit of time, that an initial state 
where two consecutive sites are occupied by molecules of type (a + £, [3 + m) end 
up occupied by molecules of type (a,/3). The rate u>o,o(«,/3) gives the probability, 
per unit of time, that the configuration of two neighbouring sites (a,/?) are changed. 
From these definitions the conservation of probabilities implies 




fc=l L 

N-l 



= - W 0fi{f3k, Pk+l)P(Pl, ■■■,pL,t) 



+ E' *>i,m(Pk, /3fc+i)P(/3i, • • • , \Pk + £}n, \Pk+i + m] N , ...,/3 L ;t) (10) 



£,m=0 




(11) 



where r = s = is again excluded. 
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The master equation ( |I0| ) can be interpreted as a Schrddinger equation 





dt 



\V) = -H \V) (12) 



in Euclidean time, if we identify the probability P({/3},t) as the wave function. The 
Hamiltonian in ( |12|) 

H=J2H J (13) 

3=1 

acts in a Hilbert space of dimension N L , while Hj acts in the subspace £g> V^ +1 ^ 
and is given by 

//, I) 'Ij. (14) 

where 

N-l N-l 

T j = E' w e,m(aif3)E a ' [a+e]N ®E^ +m]N (15) 

f,m=0 a,/3=0 
N-l 

Uj = E w 0jQ {a,(3)E aa ®E^. (16) 

o,/3=0 

The above Hamiltonian describes the time evolution from an initial probability ^>o({P},t 
0) with dynamical processes given by the transition rates w^ m {a,(3). 

Here are some examples of processes in which the HjS of Eq. ([T3|) coincide with 
the generators ej of the Hecke algebra H n (q). 

a) diffusion processes: we consider (N — 1) types of molecules (Ax, . . . , Ajv-i) 
which diffuse to the right: 

A b + (fr^(& + A b (6 = 1,...,JV-1) (17) 
with equal rates Tr and to the left: 

(fr + A b ^A b + (& (6=1,...,JV-1) (18) 
with equal rates Tl- Choosing the time scale in such a way that the diffusion constant 



D = ^T L T R = 1 and introducing q = yT R /T L} up to an equivalence transformation, 
we obtain Hj = ej corresponding to the quotient (2, 0). 

b) diffusion and interchange processes: to the diffusion processes (|T7| ) and 
( |T%D we add the interchange to the right processes: 

A b + A c -> A c + A b (b^c,b> c) (19) 

with rates T R and with rates Ti if b < c. In this case the quotient is (N, 0) and 
the obtained representation of H n (q) is hermitian (the same applies to the diffusion 
processes of case a) . 
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c) diffusion, interchange and coagulation processes: we now add to the 
processes (p~7| - p79D only one reaction for the molecule A\ (the molecules are ordered!) 



A 1 + A 1 -> + A 1 (20) 

with a rate and the processes where and A\ interchange their position in the final 
state with rate Y L . In this case the quotient is (iV — 1, 1) and the representation of H n 
is non-hermit ian. As a rule, non-hermitian representations occur for quotients (P, M) 
with M/0 i.e. when the representative Hamiltonians if ( p > A/ ) are supersymmetric. 
One obtains new quotients of H n generalizing the reaction (p0|) by taking 



A r + A r -> A s + A r (rate T R ) (21) 

A r > + A r > -> A s > + A r > (rate T R ) (22) 

where s — r ± 1, s' — r' ± 1, and r 7^ r', s'; s 7^ r', s'. In this way one obtains the 
quotient (N — 2,2). Taking more reactions like in Eq.([H],^) one can get arbitrary 
quotients (P, M). 

d) diffusion-polymerisation (modulo N) processes: in this case besides the 
diffusion (Eqs. [1T| - [T8D , one considers the processes: 

A r + A s -> + A [r+s]N (rate T R , [r + s] N ^ 0) (23) 

A r + A s -> A [r+s]jv + (rate T L , [r + s]^ ^ 0) (24) 

A r + A s ^0 + (rate Y R + T L , [r + s] N = 0) (r, s = 1, 2, . . . , N - 1). (25) 

For all these processes one obtains the quotients (1,1). Let us remark that the 
processes described by Eq. (|3|-|25|) were already considered in Ref. |§ and called 



"P-poly" for irreversible diffusion cluster-cluster aggregation processes. Using prob- 
abilistic methods several exact results were obtained (see also Ref.|10[) concerning 
the concentrations of the different types of molecules. We believe that the profound 
reason which explains why these results could have been obtained is that the (1, 1) 
quotient has a spectrum which can be obtained using free fermion methods [TI|) (see 
also @). 

e) diffusion, interchange and annihilation processes: because of interchange 
processes, the labels of the molecules are ordered; we can now add the following two 
processes: 

A m + A m -> A m+1 + A m+1 (rate Ai) (26) 

A m + A rn -> A m _ x + A m „! (rateA_ x ), (27) 

with Ai + A_i = + T^. In this way we obtain again the quotient (2, 1). 

The reader might wonder how the above H n (q) quotients were obtained. From the 
work done in Ref.0 we have developped a certain intuition for chemical processes 
up to N = 3. We have then used the computer to check that no other solutions 
exist for N = 2 and 3 and generalized the solutions obtained for A up to 3 to all 
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N and checked them. There are certainly other solutions which one can discover 
studying systems with N = 4 and certainly examples of the quotients (2, 2) different 
from the the one given in Eq . ( ^I|p2l) will show up. We hope that the present effort 
to connect representatives of H n (q) and chemical processes will help to compute 
various correlation functions for the latter. One immediate consequence of the above 
mentioned connecxion is that all processes with q ^ 1 are massive (the concentrations 
for example have an exponential fall-off in time) and all those with q = 1 are massless 
(algebraic fall-off). 
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